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ABSTRACT 


-^The  forces  acting  on  yawed  smooth  and  rough  circular 
cylinders  in  oscillating  flow  have  been  investigated  for 
the  purpose  of  determining  the  appropriate  force-transfer 
coefficients  and  the  applicability  of  the  ^'independence 
principled  The  results  have  shown  that  the  flow  about 
each  cylinder  is  unique  and  the  independence  principle 
does  not  hold  true.  In  addition,  a  comprehensive  sensi¬ 
tivity  analysis  has  been  carried  out  to  investigate  the 
effects  of  undetected  environmental  disturbances  on  the 


force- transfer  coefficients  and.  the  limits  of  probable 
errors  have  been  delineated.. 
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I.  INTRODUCTION 


The  study  of  time  dependent  flow  about  yawed  cylinders  is 
a  topic  which  presently  has  great  interest,  both  theoretically 
and  practically.  Of  the  great  variety  of  time-dependent  flows 
about  yawed  cylinders,  the  wave  motion  and  sinusoidally  oscil¬ 
lating  flow  are  of  such  importance  as  to  warrant  immediate 
investigation.  For  lack  of  adequate  information  concerning 
time-dependent  flows,  industry  is  often  forced  to  adopt 
steady  flow  relations  to  time-dependent  flow  situations. 

This  research  proposes  to  investigate  the  forces  acting  upon 
yawed  circular  cylinders  immersed  in  sinusoidally  oscillat¬ 
ing  uniform  flow  and  to  examine  the  validity  or  limitations 
of  the  present  methods  of  analysis. 

In  1950,  Morison,  et  al.  [1]  introduced  an  equation  for 
calculating  the  in-line  force  due  to  unbroken  surface  waves 
acting  on  a  vertical  pile.  For  a  cylinder  of  diameter  D, 
the  force  per  unit  length  is  expressed  as: 

F  -  J  0  D  Cd  U|U[  +  0.25  p  .D2  Cm  jjj£  (1) 

where  U  represents  the  incident  flow  velocity?  C^,  the  drag 

coefficient;  and  C  ,  the  inertia  coefficient.  The  coeffi- 

m 

cients  C  and  Cj  were  considered  to  be  time-invariant  and 
m  a 

constant  along  the  length  of  the  cylinder.  Morison 's  paper 


was  intended  as  a  preliminary  report  with  follow-on  studies 
to  be  done  on  other  structures  in  various  wave  actions. 
Morison,  et  al.  did  not  consider  the  contributions  of 
transverse  forces  and  vortex  shedding  in  the  calculation 
of  the  in-line  forces. 

Attempts  have  been  made  to  extend  this  rather  simple 
relation  to  far  more  complex  situations  such  as  combined  wave 
and  current  flow ,  hydro-elastic  oscillations,  and  yawed 
cylinders.  Heideman,  et  al*  [2]  studied  Morison' s  equation 
using  ocean  data  and  concluded  that  Morison' s  equation  was 
satisfactory  for  normal  cylinders.  The  utilization  of  an 
equation  unproven  even  under  ideal  flow  situations  to 
correlate  ocean  data  so  as  to  prove  the  validity  of  the 
equation  is  not  a  very  meaningful  exercise.  Sarpkaya  and 
Isaacson  [3]  reviewed  in  detail  the  numerous  studies  regard¬ 
ing  the  application  of  the  Morison  equation. 

It  is  obvious  that  there  is  a  great  need  for  an  ideali¬ 
zation  of  the  problem,  or  an  experiment  which  is  more 
manageable  (e.g.,  a  sinusoidally  oscillating  uniform  flow). 
Only  in  this  way  can  all  the  complex  interactions  be 
separately  taken  into  account. 

Engineers  faced  with  the  problem  of  dealing  with  wave 
forces  on  yawed  members,  having  no  other  recourse,  drew 
upon  previous  work  with  steady  flow.  Hoerner  [4]  proposed 
the  "independence  principle",  which  stated  that  the  normal 
pressure  forces  are  independent  of  the  tangential  velocity 


for  subcritical  values  of  Ren,  where  Ren  is  the  Reynolds 
number  based  upon  the  component  of  the  flow  velocity  normal 
to  the  cylinder.  This  principle  allowed  the  neglecting  of 
the  tangential  components.  Bursnall  and  Loftin  [5]  found 
that  the  independence  principle  does  not  apply  to  the  critical 
and  transcritical  flow  regimes.  Norton,  et  al.  [6]  found  that 
the  independence  principle  does  apply  to  the  postcritical  as 
well  as  subcritical  flow,  but  not  to  the  critical  and  trans¬ 
critical  regions  in  between.  Thus,  recent  research  has  shown 
that  the  independence  principle  applies  when  the  boundary 
layer  is  wholly  laminar  (Hoemer)  or  wholly  turbulent  (Norton)  , 
but  its  use  in  the  critical  and  transcritical  regions  is 
uncertain. 

The  designers  of  offshore  structures  were  thus  led  to 
adopt  the  independence  principle  for  the  wave  force  calcula¬ 
tions  and,  in  doing  so,  to  generalize  the  Ilorison  equation. 

This  has  been  accepted  practice  in  industry  for  lack  of  a 
better  relation.  It  can  be  asked  what  would  be  anticipated 
for  oscillating  flow  or  waves  on  the  basis  of  what  is  known 
for  steady  flow.  The  instantaneous  Reynolds  number  in  such 
a  flow  will  vary  from  -Remsiv  to  +Rem=v  during  a  complete  flow 
cycle.  It  could  be  postulated  that  the  boundary  layer  would, 
at  times,  be  fully  laminar;  at  other  times  fully  turbulent; 
and  the  rest  of  the  time  be  in  transition.  In  light  of  this, 
it  is  rather  doubtful  that  the  independence  principle  applies 
at  all  to  oscillating  or  wave  flow. 
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The  study  of  forces  acting  on  yawed  cylinders  can  be 
investigated  either  by  oscillating  a  yawed  cylinder  in  a 
tank,  using  small  amplitude  waves  in  a  laboratory  channel, 
oscillating  the  flow  about  a  fixed  cylinder,  or  by  using 
ocean  data  if  available.  Oscillating  the  cylinder  has 
proven  to  be  impractical  because  of  difficulties  in  account¬ 
ing  for  the  inertial  force  acting  on  the  body,  producing 
repeatable  oscillations,  controlling  spurious  vibrations  in 
the  system,  and  in  measuring  the  in-line  and  traverse  forces 
simultaneously .  Waves  have  relatively  more-complex  flow 
kinematics  due  to  the  orbital  motion  of  particles  and  the 
decay  of  wave  amplitude  with  depth.  Ocean  data  are  not 
available  for  such  a  study.  Oscillation  of  the  flow  past 
the  cylinder  has  proved  to  be  the  most  practicable  method. 

By  using  a  sinusoidally  oscillating  uniform  flow  (which  can 
be  fairly  easily  produced)  the  kinematics  of  the  flow  are 
accurately  known.  Sarpkaya  and  Isaacson  [3]  discuss  the 
similarities  and  differences  between  such  a  flow  and  the 
flows  described  by  wave  theory. 

The  primary  objective  of  this  investigation  was  to  study 
the  forces  exerted  by  a  sinusoidally  oscillating  uniform  flow 
on  yawed  circular  cylinders  to  determine  whether  or  not  the 
independence  principle  is  applicable.  If  the  principle  holds, 
the  force  transfer  coefficients  calculated  by  Fourier  analysis 
using  the  normal  velocity  component  should  reduce  identically 
to  the  normal  cylinder  case  at  corresponding  values  of  K,  Re, 
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and  k/D.  If  the  independence  principle  does  not  apply,  it 
is  desired  to  determine  what  the  coefficients  are  as  functions 
of  yaw  angle,  roughness  ratio  k/D,  Re,  and  K.  It  would  also 
be  necessary  to  determine  how  well  Morison’s  equation  works 
with  the  new  coefficients. 

With  the  foregoing  objectives  in  mind  a  detailed  investi¬ 
gation  was  undertaken  using  smooth  and  rough  cylinders  of 
nominal  6-inch,  4^— inch,  and  3-inch  diameters  at  yaw  angles 
of  45  degrees,  60  degrees,  and  90  degrees  (as  measured  between 
the  ambient  flow  direction  and  the  cylinder  axis) .  Part  of 
this  effort  has  already  been  reported  [7],  The  present  study 
extends  the  previous  work,  accomplishes  the  original  objec¬ 
tives  of  the  investigation,  and  takes  a  much  closer  look  at 
the  sensitivity  of  the  modified  Morison  equation  and  of  the 
measured  forces  to  random  external  disturbances  encountered 


in  nature. 


EXPERIMENTAL  APPARATUS  AND  CYLINDERS 


The  experiments  were  carried  out  in  the  large  oscillating- 
flow  tunnel  at  the  Naval  Postgraduate  School.  The  tunnel  was 
basically  unchanged  from  the  configuration  described  in  [8 
and  9].  Mounting  pads  for  the  strain  gage  housings  were 
placed  on  opposite  sides  of  the  tunnel  in  such  a  manner  so 
as  to  accommodate  the  yawed  cylinders  (Fig.  1) .  The  mounting 
pad  on  one  side  of  the  tunnel  was  fitted  with  a  slide  so  that 
the  pad  could  be  adjusted  small  distances  to  exactly  match 
the  cylinder  length.  The  force  transducers  and  housings  were 
unchanged  from  those  described  in  [9]. 

The  cylinders  also  were  similar  to  those  previously 
reported,  except  that  the  ends  were  cut  exactly  at  the 
desired  angle  of  yaw.  The  cylinder  ends  were  parallel  to 
the  tunnel  walls  with  a  1/32-inch  gap  at  each  end.  This  gap 
was  filled  with  a  soft  foamy  material  glued  to  the  ends  of 
the  cylinder.  Ball  bearings  were  mounted  in  the  ends  of  the 
cylinder  with  the  outer  bearing  faces  flush  with  the  face  of 
the  cylinder. 

The  cylinders  mounted  in  the  tunnel  were  supported  by 
the  pins  in  the  force  transducers.  The  pins  were  perpendicu¬ 
lar  to  the  tunnel  wall.  Rubber  washers  were  placed  behind 
the  pins  so  that  the  lateral  motion  of  the  cylinders  (normal 
to  the  tunnel  wall)  was  restrained,  but  without  compressing 
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the  cylinder.  The  cylinders  were  free  to  rotate,  but  a  small 
clip  was  placed  on  one  end  of  the  cylinder  to  restrict  the 
rotation  to  a  very  small  amount  so  that  the  yawed  end  of  the 
cylinder  did  not  strike  the  tunnel  wall  as  it  rotated.  Thus, 
as  mounted,  the  cylinder  was  laterally  constrained,  but  free 
to  rotate  within  a  small  range.  Careful  calibration  and 
extensive  testing  proved  that  this  mounting  system  allowed 
accurate  and  repeatable  recording  of  the  normal  and  trans¬ 
verse  forces  acting  on  the  cylinder. 

Calibration  was  first  conducted  in  a  vertical  direction, 
since  yaw  would  not  affect  the  total  force  in  this  direction. 
Weights  were  hung  from  the  center  of  the  horizontal  cylinder, 
and  calibration  factors  were  obtained  by  converting  the  sum 
of  the  signals  of  two  force  transducers  to  pounds/mm  and 
pounds/volt.  The  strain  gages  were  found  to  be  linear 
throughout  the  range  of  expected  forces.  The  horizontal 
forces  normal  to  the  cylinder  are  related  to  the  streamwise 
force  by  the  sine  of  the  yaw  angle  a  (as  measured  between  the 
cylinder  axis  and  ambient  flow) .  a  system  of  supports  and 
pulleys  allowed  applying  horizontal  force  at  the  center  of 
the  cylinder  and  normal  to  it  (Fig.  1) .  The  sum  of  the 
outputs  of  the  two  transducers  was  recorded  on  the  strip 
chart  recorder,  establishing  the  necessary  relationship 
between  the  normal  in-line  force*  F  and  the  total  electrical 

#The  force  in  the  plane  defined  by  the  velocity  vector  and 
the  cylinder  axis. 
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output  of  the  gages.  The  same  gages  were  used  to  measure 
the  lift  force  or  the  in-line  force  by  rotating  the  gages 
90  degrees.  It  is  easy  to  show  that  the  calibration  factor 
of  the  lift  force  (here  in  the  vertical  direction,  up  or 
down) ,  expressed  in  terms  of  pounds  or  millimeter  deflection, 
is  sina  times  the  calibration  factor  for  the  in-line  force. 
This  proved  to  be  true  experimentally  also  and  demonstrated 
independently  the  validity  of  the  measurement  technique. 

Raines  [7]  carried  out  experiments  with  the  following 
cylinders : 

90  degrees,  5.905  inch,  rough  and  smooth  (nominal 
size  6  inches) 

60  degrees,  6.005  inch,  rough 

45  degrees,  5.923  inch,  rough  and  smooth 

45  degrees,  4.42  inch,  rough  and  smooth  (nominal 
size  4.5  inches) 

45  degrees,  3.00  inch,  rough  and  smooth  (nominal 
size  3  inches) 

The  present  experiments  were  carried  out  with  the 
following  cylinders : 

60  degrees,  5.965  inch,  smooth  (nominal  size 
6  inches) 

60  degrees,  4.450  inch,  smooth  and  rough 
(nominal  size  4.5  inches) 

60  degrees,  2.985  inch,  smooth  and  rough 
(nominal  size  3  inches) 

Relative  sand  roughness  of  k/D  =  1/100  was  used  for  all 
rough  cylinders. 
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In  addition,  the  foregoing  data  were  complemented  with 
data  reported  in  110].  These  data,  taken  using  3  inch  and 
4.5  inch  nominal  diameter  cylinders  at  90  degree  yaw  angle, 
were  used  in  the  analysis  and  comparison  of  the  entire  data. 
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III.  DATA  ACQUISITION  AND  PROCESSING 


Data  acquisition  was  accomplished  using  an  HP3052A 
automatic  data  acquisition  system.  The  system  was  configured 
as  described  in  [9] . 

The  data  acquisition  system  consisted  of  an  HP3495A 
20  channel  scanner  and  an  HP3437A  system  voltmeter,  both 
controlled  by  an  HP9845B  desk  top  computer.  The  in-line 
or  transverse  signal,  originating  as  a  voltage  from  the 
strain  gages,  was  amplified  and  sent  to  the  strip  chart 
recorder.  This  amplified  signal  was  also  sent  to  the  scan¬ 
ner  as  one  channel  of  data.  The  flow  amplitude  signal 
originated  from  a  differential  pressure  transducer,  and  was 
similarly  amplified  and  sent  to  both  the  strip  chart  recorder 
and  the  scanner  as  another  channel  of  data.  Strip  chart  data 
was  maintained  as  a  visual  record  of  the  flow  amplitude  and 
force.  The  scanner  switched  between  the  two  channels  at 
designated  intervals  when  triggered  by  the  voltmeter.  The 
voltmeter  read  each  one  of  these  voltage  values  and  trans¬ 
ferred  them  to  the  computer  for  immediate  calculation  of 
the  governing  parameters. 

At  least  six  cycles  of  in-line  force  data,  recorded 
aoproximately  70  times  per  cycle,  were  acquired  with  the 
system.  The  data  were  then  averaged  to  give  one  cycle  from 
which  the  governing  parameters  were  calculated.  Ten  cycles 
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of  transverse  force  data,  recorded  at  the  same  rate,  were 
acquired  with  the  system.  The  transverse  data  was  not 
averaged  prior  to  calculating  the  governing  parameters. 

The  use  of  the  HP3052A  system  proved  highly  beneficial 
for  a  number  of  reasons.  The  calculated  parameters  were 
available  on  a  real  time  basis,  allowing  for  detailed  study 
of  anomalies  in  the  force  patterns.  Taking  data  over  a 
large  number  of  cycles  removes  more  of  the  randomness  of 
the  flow  forces.  This  was  particularly  valuable  for 
analyzing  the  transverse  force  data. 
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IV.  GOVERNING  PARAMETERS 

Data  reduction  for  the  in-line  forces  is  based  on 
Morison's  equation  using  the  normal  component  of  velocity. 
The  derivation  of  the  governing  parameters  is  described  in 
detail  in  [3  and  9] .  Here  only  a  brief  description  of  the 
parameters  will  be  given. 

For  a  sinusoidally  varying  flow,  Eq.  (1)  may  be 
written: 
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where  the  flow  is  represented  by  U  =  -Umcos9  with  6  =  2irt/T 
The  Fourier  averages  of  Cm  and  Cd  are  given  as  [1] : 
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where  Fm  represents  the  measured  force. 
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As  discussed  in  [9] ,  it  would  be  beneficial  to  remove 
the  cycle-to-cycle  variations  in  the  force.  This  was  accom¬ 
plished  by  using  the  normalized  RMS  value  of  the  actual  force 
(C  > .  Two  new  coefficients  were  defined  from  this: 


c.  =  c./c 

dr  d/  rms 


C  ,  =  C  ./C 

dr  d  amrs 


(5) 


and 


c  =  c  /c  ,  c  =  c  /c 

mr  nr  rms  mr  m  arms 

Using  Eq.  (2) ,  the  rms  value  of  the  calculated  force 
can  be  shown  to  be : 
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Transverse  force  data  were  represented  by  the  normalized 
rms  value  of  the  measured  lift  force  •  Previously, 

some  of  the  lift  data  were  analyzed  [7]  using  a  Fourier 
transform  to  obtain  frequency  information.  In  the  present 
study,  this  effort  was  discontinued  because  the  frequency 
information  so  far  obtained  was  more  than  sufficient  to 


draw  general  conclusions  regarding  the  frequency  content  of 
the  lift  force. 


V.  EXPERIMENTAL  RESULTS 

A.  PRESENTATION  OF  THE  RESULTS 

The  force-transfer  coefficients  (C  ,  C.,  CaT__,  C _ ,  C-  , 

m  a  arms  mr  dr 

C-,  rTng)  are  presented  in  graphical  form  in  terms  of  the 
governing  parameters  K,  S,  a,  and  k/D  in  Figs.  2  through  7. 

The  new  as  well  as  the  previously-obtained  data  [7  and 
10]  have  been  cross-plotted  in  Figs.  8  through  13  in  order 
to  exhibit  the  effects  of  the  governing  parameters. 

B.  DISCUSSION  OF  THE  RESULTS 

1.  Inertia  and  Drag  Coefficients 

The  first  indication  of  the  fact  that  the  independence 
principle  is  not  applicable  to  yawed  cylinders  in  oscillating 
flow  is  seen  in  Fig.  8a  where  Cm  is  plotted  as  a  function  of 
K  for  8  =  4040,  k/D  =  0.00,  and  a  =  90 ,  60,  and  45  degrees. 
Evidently,  had  the  independence  principle  been  valid,  the 
three  sets  of  data  corresponding  to  three  yaw  angles  would 
have  collapsed  into  a  single  set,  thus  rendering  the  Cm 
versus  K  relationship  independent  of  the  yaw  angle  (at  least 
within  the  range  of  the  test  parameters  and  experimental  errors) . 
This  conclusion  is  more  universal  than  it  may  appear  at  first 
sight  on  the  basis  of  Fig.  8a.  In  fact.  Figs.  8b  through  8f, 
presenting  Cm  versus  K  for  both  smooth  and  rough  cylinders  of 
6,  4.5,  and  3  inch  nominal  diameters,  show  that  C  varies  with 
the  yaw  angle  as  it  did  in  Fig.  8a.  Thus,  one  may  conclude 
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that  the  independence  principle  is  not  valid  at  least  for  the 
inertia  coefficient.  The  recognition  of  this  fact  is  impor¬ 
tant  for  two  reasons.  Firstly,  the  inertial  force  is  only 
part  of  the  total  force.  Secondly,  the  significance  of  the 
inertial  force  depends  on  K.  Thus,  the  failure  of  the 
inertia  coefficient  to  obey  the  independence  principle  may 
not  by  itself  be  very  important  if  the  inertial  force  is 
small  relative  to  the  total  force  for  the  particular  case 
under  consideration.  Consequently,  it  is  important  to 
examine  the  variation  of  the  drag  coefficient  and  its 
relation  to  the  total  force  and  the  independence  principle. 

Figure  9b  shows  the  variation  of  Cd  as  a  function  of 
K  for  3  =  2345,  k/D  =  0.00,  and  ot  =  90,  60,  and  45  degrees. 
Clearly,  the  variation  of  Cd  with  K  and  the  yaw  angle  is  far 
more  complex  than  that  of  C^.  Figures  9a  through  9c, 
depicting  C^  versus  K  for  the  three  smooth  cylinders,  show 
that  the  independence  principle  does  not  hold  true  except  at 
relatively  large  values  of  K.  The  implication  of  this  result 
is  that  when  the  vortex  motion  about  the  body  becomes  ex¬ 
tremely  turbulent  and  incoherent,  as  it  would  at  high  K 
values,  the  overall  effect  of  yaw  becomes  less  significant. 

It  is  also  in  this  region  that  the  force  is  drag  dominated 
and  the  drag  force  may  be  calculated  in  terms  of  the  force 
acting  on  the  projected  area.  Though  the  inertia  coefficient 
exhibits  strong  variation  with  yaw  even  at  high  K  values, 
this  variation  is  of  academic  interest  only  since  the 
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contribution  of  the  inertial  force  to  the  total  force  is 
rather  insignificant. 

Roughness  plays  an  important  and,  oddly  enough,  a 
unifying  role  on  the  drag  coefficient.  This  is  somewhat 
anticipated  on  the  basis  of  the  fact  that  roughness  precip¬ 
itates  earlier  transition  and  greater  incoherence  in  the 
flow  about  a  yawed  cylinder,  thereby  creating  conditions 
more  favorable  to  the  independence  principle.  Figures  9d 
through  9f  show  that  the  independence  principle  for  rough 
cylinders  is  almost  valid  with  the  exception  of  the  drag- 
inertia  dominated  regime  (8<K<20). 

It  is  clear  from  the  foregoing  discussion  of  the 
inertia  and  drag  coefficients  that  the  normal  force  acting 
on  a  yawed  cylinder  is  significantly  underestimated  through 
the  use  of  Morison's  equation,  independence  principle,  and 
the  drag  and  inertia  coefficients  appropriate  to  the  normal 
cylinder.  Only  a  narrow  range  of  K  values  (from  about  9  to 
17),  in  the  drag- inertia  dominated  regime,  that  the  drag 
coefficient  for  a  normal  cylinder  exceeds  that  for  a  yawed 
cylinder.  However,  the  large  increase  in  the  inertia 
coefficient  of  the  yawed  cylinder  in  the  same  region  of  K 
values  more  than  offsets  the  relatively  small  decrease  in 
its  drag  coefficient. 

The  absence  of  the  inertia  crisis,  the  large  increase 
in  the  inertia  coefficient,  the  relative  decrease  in  the  drag 
and  lift  coefficients,  all  in  the  drag- inertia  dominated 
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regime,  of  yawed  cylinders  are  a  consequence  of  the  decrease 
of  the  spanwise  coherence  of  vortices.  In  general,  the 
smaller  the  spanwise  coherence,  the  larger  the  Cm  in  this 
region.  It  appears  that  the  independence  principle  or  the 
cosine  law,  as  it  is  sometimes  called,  is  a  gross  simplifica¬ 
tion  of  the  behavior  of  flow  in  the  near  wake.  One  may, 
therefore,  conclude  that  the  Fourier-averaged  drag  and  inertia 
coefficients,  based  on  Morison’s  equation,  are  unique  for 
each  angle  of  yaw,  Reynolds  number,  Keulegan-Carpenter  number, 
and  the  relative  roughness. 

2.  RMS  Values  of  the  Inertia  and  Drag  Coefficients 

It  has  already  been  shown  [7]  that  the  normalized  rms 
values  of  the  inertia  or  drag  coefficients  remove  signifi¬ 
cantly  the  scatter  from  the  data  and  establish  clearer 
relationships  between  these  coefficients  and  the  governing 
parameters.  The  same  rms  coefficients  are  used  here  to 
examine  the  effect  of  yaw  and  experimental  error.  Figures 

11a  through  Ilf  show  C _  =  C  /C  _  as  a  function  of  K. 

Clearly,  Cmr  is  independent  of  Reynolds  number  (or  3)  ,  yaw 
angle,  and  the  roughness  for  K  smaller  than  about  8.  For 
smooth  cylinders ,  Cmr  is  nearly  independent  of  yaw  for  K 
larger  than  about  20.  Additionally,  for  rough  cylinders  Cmr 
is  least  dependent  on  yaw  for  relatively  low  values  of  Re. 

The  drag  coefficient  C^r  exhibits  similar  character¬ 
istics  (see  Figs.  12a  through  12f) .  The  scatter  in  Cdr  is 
considerably  less  than  that  in  Cd  primarily  because  of  the 


31 


fact  that  the  variations  in  Cd  as  a  consequence  of  the  random 

disturbances  imposed  on  the  measured  force  are  smoothened  out 

by  the  variations  in  the  corresponding  Carmg  values.  This 

fact  has  been  substantiated  through  numerous  examples  by 

examining  the  relative  variations  of  C.,  C  ,  and  C  .  As 

am  arms 

K  increases  Cdr  reaches  its  theoretical  value  of  1.63,  as 
seen  from  Eq.  (7) .  Finally,  roughness  tends  to  decrease  the 
dependency  of  Cdr  on  Re  (see  Figs.  12d  through  12f) . 

3.  Transverse  (Lift)  Force 

The  lift  force  data  are  presented  in  terms  of  the 
normalized  rms  values  of  the  measured  lift  force  in  Figs.  13a 
through  13c  for  the  three  smooth  cylinders.  These  figures 
show  that  C-,  T^ng  is  largest  for  the  90-degree  cylinders,  as 
one  would  expect,  and  decreases  rapidly  with  yaw.  However, 
the  variation  of  Clrmg  with  the  yaw  angle  (60  and  45  degrees) 
is  not  as  large  as  that  with  the  Reynolds  number.  The  lift 
coefficient  increases  rapidly  as  K  increases  from  about  6 
to  12.  Then,  ^  again  decreases  rapidly  due  to  the 
increased  lack  of  coherence  of  the  vortices.  Similar 
observations  may  be  made  for  rough  cylinders  (see  Figs.  13d 
through  13f) .  It  suffices  to  note  that  the  yaw  has  an 
immediate  effect  on  the  coherence  of  vortices  and  on  the  lift 
coefficient,  particularly  in  the  drag- inertia  dominated 
regime.  For  large  values  of  K  or  Re,  the  irregularity  of 
the  vortex  shedding  even  for  a  90-degree  cylinder  obscures 
the  yaw  effect.  It  is  then  rather  difficult  to  distinguish 
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between  the  lift  coefficients  of  cylinders  at  various  yaw 
angles.  In  the  ocean  environment,  one  is  more  likely  to  find 
a  confused  vortex  shedding  with  very  small  coherence  length. 
Thus,  one  may  safely  state  that  the  inclined  members  subjected 
to  the  wave  motion  in  the  ocean  environment  will  experience 
considerably  less  lift  force  than  those  normal  to  the  waves. 
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VI.  SENSITIVITY  ANALYSIS 


A.  PHASE  SHIFT  EFFECTS 

It  has  been  known  since  the  inception  of  Morison’s 
equation  in  1950  that  the  phase  change  between  the  maximum 
force  and  the  maximum  velocity  may  lead  to  varying  degrees 
of  error,  depending  on  the  governing  parameters.  This  phase 
shift  may  arise  due  to  electronic  or  numerical  filtering  of 
the  force  and  velocity  data  or  due  to  variations  in  the 
coherence  of  the  vortices  along  a  given  separation  line. 

Even  though  it  is  not  yet  possible  to  attribute  a  given 
fraction  of  the  estimated  error  to  a  specific  cause,  it  is 
most  desirable  to  examine  in  detail  the  isolated  effects  of 
such  things  as  an  imposed  phase  shift  and/or  a  random 
disturbance  of  a  given  amplitude  and  frequency.  Eventually, 
such  information  will  show  the  region  of  sensitivity  of  the 
force  coefficients  and  Morison’s  equation  to  the  disturbances 
imposed  by  the  environmental  conditions  and  will  lead  to  a 
more  rational  interpretation  of  the  causes  of  large  scatter 
observed  in  the  results  of  ocean  tests. 

The  effects  of  phase  shift  on  the  force-transfer 
coefficients  have  been  examined  as  follows.  For  a  given 
K,  0,  k/D,  and  yaw  angle  the  drag  and  inertia  coefficients 
have  been  calculated  from  Eqs.  (3)  and  (4) .  These 
coefficients  correspond  to  cases  where  there  was  no  phase 


shift  as  ascertained  by  the  manner  in  which  the  experiments 
were  carried  out.  Then  the  force  was  recalculated  through 
the  use  of  Morison's  equation  and  the  above  coefficients. 

The  calculated  force  was  then  given  a  phase  shift  and  a  new 
set  of  force  coefficients  were  calculated  through  the  use  of 
Eqs.  (3)  and  (4).  The  method  was  repeated  for  representa¬ 
tive  values  of  K  and  reasonable  phase  angles. 

Figures  14a  through  14e  show  the  results  for  five 
coefficients.  Clearly,  the  phase  angle  has  significant 
effect  on  all  of  these  coefficients  with  the  exception  of 
Carms  Concentrating  attention  on  Cm  and  Cd  it  is  noted 
that  Cra  decreases  and  increases  with  increasing  phase 
angle  (where  the  force  curve  is  shifted  forward  in  time) . 

The  relative  rate  of  decrease  or  increase  of  these  two 
coefficients  certainly  depends  on  the  particular  values  of 
the  governing  parameters.  For  example,  for  K  =  6.66  (Fig. 
14a)  Cm  decreases  rather  gradually  whereas  increases 
rapidly  with  increasing  phase  angle.  The  reason  for  this 
is  that  for  values  of  K  less  than  about  8  the  fluid  force 
is  inertia  dominated  and  the  inertia  coefficient  is  not 
very  sensitive  to  small  disturbances  imposed  on  it.  In 
mathematical  terms  the  foregoing  simply  means  that  the  force 
curve  is  like  a  sine  curve  and  that  when  it  is  shifted  right 
or  left  by  some  amount  and  then  multiplied  with  sin9,  as  in 
Eq.  (3) ,  the  integral  could  not  be  too  sensitive  to  the 
phase  difference.  The  drag  coefficient,  on  the  other  hand, 
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is  extremely  sensitive  to  small  phase  shifts  in  the  inertia 
dominated  region  simply  because  a  small  phase  3hift  in  a 
"sine-like”  curve  multiplied  with  Icos9|cos0,  as  in  Eq.  (4), 
produces  large  changes. 

In  the  drag-inertia  dominated  regime  (K  ranging  from 
about  8  to  20)  both  and  Cm  are  affected,  as  seen  in  Figs. 
14c  and  14d.  In  the  drag  dominated  regime,  however,  Cd  is 
hardly  affected  by  phase  shifts;  whereas  Cm  shows  a  dramatic 
drop  with  increasing  phase  angle  (Fig.  14e) .  The  reason  for 
this  is  easily  explainable  in  a  manner  similar  to  that 
previously  presented  in  connection  with  Fig.  14a. 

The  overall  conclusion  to  be  drawn  from  the  foregoing 
is  that  there  is  no  flow  regime  where  one  or  the  other 
coefficients  is  not  effected  by  the  phase  shift.  In  fact, 
if  a  constant  phase  shift  of  ±5  degrees  were  to  be  randomly 
imposed  on  a  series  of  data  obtained  with  K  values  ranging 
from  5  to  100  one  would  have  observed  large  scatter  in  Cd 
and  small  scatter  in  for  K  smaller  than  about  8;  large 

scatter  in  both  CA  and  C„  for  K  from  about  8  to  20 ;  and 
large  scatter  in  Cm  and  small  scatter  in  Cd  for  K  larger 
than  about  20.  These  scatter  patterns  are  invariably 
observed  in  the  data  obtained  in  the  ocean  environment. 

As  noted  earlier,  not  all  of  the  phase  shift  is  due  to 
electronic  filtering  of  the  data.  Some  of  the  phase  shift 
comes  from  the  reduction  of  the  coherence  length  of  the 
vortices  by  the  disturbances  ever-present  in  the  wave  field. 
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Thus  the  purpose  of  the  foregoing  is  not  to  suggest  that 
all  phase  shifts  could  be  avoided  or  eliminated  but  rather 
to  quantify  the  effects  of  the  phase  shift  in  terms  of  the 
expected  variation  of  the  force  coefficients  in  order  to  have 
a  better  grasp  of  the  meaning  of  the  observed  scatter. 

B.  RANDOM  SQUARE  WAVE  EFFECTS 

The  experiments  conducted  under  controlled  laboratory 
conditions  are  immune  to  practically  all  of  the  disturbances 
encountered  in  the  ocean  environment.  Consequently,  the 
vortices  have  larger  coherence  lengths  and  these  give  rise 
to  maximum  values  of  the  semi- peak -to- peak  force  coefficient 
Cp(spp).*  In  the  ocean  environment,  however,  the  variable 
wave  direction,  change  of  wave  velocity  with  depth,  effect 
of  secondary  waves  with  all  sorts  of  harmonics,  irregular 
surface  roughness,  and  numerous  other  stochastic  variables 
give  rise  to  varying  degrees  of  scatter  in  the  force 
coefficients  based  on  ocean  data.  It  is  important  from  the 
engineering  point  of  view  to  examine  the  collective  impact 
of  these  disturbances  on  the  force  coefficients.  It  is  under¬ 
stood  that  very  high  frequency  disturbances  imposed  on  the 
measured  force  cannot  significantly  affect  and  Cm.  On 
the  other  hand,  disturbances  which  simulate  either  in  part 

*This  does  not  imply  that  perfect  coherence  leads  to 
maximum  values  of  C3  and  Cm.  In  fact,  a  force  trace 
resulting  from  a  perfect  coherence  can  be  rendered 
incoherent  by  a  random  square  wave.  The  new  force 
trace  yields  a  larger  C  and  a  smaller  C^. 
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or  in  combination  a  "sine  like"  or  a  "cosine  like"  wave  are 


most  effective  in  causing  large  changes  in  and  C  .  This 
follows  from  the  very  definition  of  the  inertia  and  drag 
coefficients  as  in  Eqs .  (3)  and  (4).  Clearly  one  must 
examine  the  effect  of  randomly  imposed  disturbances  of 
given  frequency  and  phase  shift  on  the  variation  of  the 
force  coefficients.  Such  an  extensive  analysis  has  been 
carried  out  in  a  manner  described  as  follows.  A  frequency 
ratio  of  fgw/ft  was  chosen  where  f  is  the  frequency  of 
the  square  wave  and  ffc  is  the  frequency  of  the  periodic 
oscillation  of  the  flow  in  the  tunnel.  Then  the  square  wave 
was  assigned  a  random  amplitude  and  initial  phase  shift 
through  the  use  of  a  built-in  random  number  generator.  The 
maximum  of  the  disturbance  was  taken  to  be  the  absolute  of 
Cp(spp)  (almost  equal  to  the  absolute  maximum  of  the  measured 
force) .  This  procedure  varied  the  amplitude  of  each  disturb¬ 
ance  by  multiplying  a  random  number  between  +1  and  -1  with 
Cp(spp).  Calculations  have  been  performed  1000  times  for  a 
given  frequency  ratio  and  force  trace.  Each  calculation 
started  with  a  randomly  selected  phase  angle  and  followed 
with  square  wave  fluctuations  whose  amplitude  varied  randomly 
as  noted  above  and  whose  frequency  was  fixed  by  the  frequency 
ratio.  Figures  15a  through  15e  show  sample  plots  of  the 
original  force,  the  disturbed  force,  and  the  square  wave 
disturbance.  These  figures  are  presented  here  for  the 
purpose  of  illustration  and  not  necessarily  for  the 


i 
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quantification  of  the  boundaries  of  the  sensitivity 
envelope. 

Figure  16a  shows  the  results  obtained  in  the  manner 
described  previously.  The  two  symbols  in  this  figure 
represent  the  relative  changes  in  Cm  and  C^.  The  computer 
program  was  directed  to  identify  the  particular  cases  in 
1000  runs  which  gave  the  largest  changes  in  and  the 
corresponding  changes  in  C^.  It  is  important  to  emphasize 
that  the  change  in  in  Fig.  16a  is  not  the  maximum 
possible  change  that  could  have  taken  place  with  a  suitable 
square  wave  disturbance.  Thus,  Fig.  16a  shows  that  there 
are  certain  frequency  ratios  which  cause  a  large  change  in 
Cm  with  very  little  change  in  the  corresponding  C^.  As 
expected,  disturbances  with  a  frequency  ratio  of  1.0  yield 
the  largest  AC^/C^  since  this  nearly  corresponds  to  a  sine 
function.  Surprisingly  enough,  high  frequency  disturbances 
have  very  little  effect  on  the  inertia  coefficient,  at  least 
in  the  drag-inertia  dominated  regime  (see  Fig.  16a).  Thus, 
one  can  filter  out  such  high  frequency  disturbances  without 
materially  affecting  the  determination  of  the  force-transfer 
coefficients . 

A  similar  calculation  has  been  performed  with  the  force 
trace  corresponding  to  Fig.  16a  to  determine  the  maximum 
change  in  and  the  corresponding  change  in  Cm  for 
representative  frequency  ratios.  Figure  16b  shows  the 
results  obtained  in  this  manner.  Several  facts  are 
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immediately  apparent.  is  not  as  sensitive  as  to  the 

randomly  imposed  disturbances  in  the  drag-inertia  regime 
even  though  ACd/Cd  is  quite  appreciable.  This  is  directly 
related  to  the  stability  of  the  integral  representing  C^, 
i.e.,  | cos 6  I  cos 6  times  a  square  wave. 

Figures  17a  and  17b  show  the  effect  of  the  disturbances 
on  Cm  and  C_j  in  the  inertia  dominated  regime  and  Figs.  18a 
and  18b  illustrate  the  similar  effects  in  the  drag-inertia 
dominated  regime.  It  is  clear  from  all  of  the  figures 
pertaining  to  the  sensitivity  of  the  drag  and  inertia 
coefficients  that  it  is  not  the  dominance  of  the  drag  or 
inertia  in  a  given  regime  that  determines  the  sensitivity  of 
a  given  coefficient.  To  be  sure,  the  said  effect  is  present 
in  Figs.  16a  through  18b.  However,  the  most  significant 
elements  of  a  disturbance  are  its  frequency  and  relative 
magnitude.  Disturbances  imposed  on  the  force  at  a  frequency 
ratio  of  about  one  or  unity  will  cause  major  changes  in 
both  and  Cm  irrespective  of  the  flow  regime.  As  to  the 
amplitude  of  the  disturbances ,  they  may  depend  on  a  number 
of  parameters  such  as  K,  Re,  turbulence  in  ambient  flow, 
interference  of  other  bodies,  etc.  In  the  present  analysis 
the  maximum  amplitude  of  the  disturbances  at  any  frequency 
were  taken  to  be  equal  to  Cp(spp)  of  the  measured  force. 
Since  this  coefficient  decreases  with  increasing  K,  it  is 
easy  to  understand  that  at  high  K  values  ACm/Cm  or  AC^/C^ 
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will  result  in  values  as  large  as  those  encountered  at 
smaller  K  values  (see  Figs.  16a  through  18b). 

The  significance  of  these  results  will  be  illustrated 
with  a  numerical  example.  Assume  K  =  14.08,  Beta  =  4040 
(as  in  Fig.  16a),  for  which  Cm  =  0.61,  Cd  =  1.63  and  Cp(spp) 
1.66.  For  a  frequency  ratio  of  f  /f.  =  1.0,  one  has 
ACm/Cm  =  2.7  and  ACd/Cd  =  0.91  from  Figs.  16a  and  16b. 
Assuming  a  disturbance  amplitude  equal  to  10  percent  of 
Cp(spp) ,  one  has: 


Cd (perturbed 
value) 


(l±ACd/Cd  *  %Cp (spp) )  *  Cd(original 

value) 


=  (1±0.91  *  0.10)  *  1.63  =  J'Jg 


and 


C 


m 


(perturbed 

value) 


(l±Ac/C_  *  %C„ (spp) )  *  C_(original 
m  ra  F  m  value) 


=  (1+2.7  *  0.10)  *  0.61  = 


The  foregoing  examples  show  that  a  random  disturbance  of 
frequency  ratio  1.0,  imposed  on  a  particular  force,  may 
have  resulted  in  relatively  large  cycle  to  cycle  variations 
of  the  force  coefficients.  Here  Cd  varies  from  a  minimum 
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value  of  1.48  to  a  maximum  value  of  1.78  and  C  varies  from 

m 

a  minimum  value  of  0.45  to  a  maximum  value  of  0.77.  These 
are  relatively  large  changes  and  show  clearly  the  importance 
of  the  following  two  factors.  First,  it  is  not  possible 
to  obtain  scatter-free  data  in  the  ocean  environment  even 
for  a  single  K  and  Re  because  of  the  ever-present  random 
disturbances  emanating  from  various  sources.  Secondly,  a 
structure  cannot  be  designed  to  withstand  the  maximum  forces 
predicted  on  the  basis  of  idealized  laboratory  experiments 
because  the  maximum  forces  may  occur  only  in  a  small  fraction 
of  the  lifespan  of  the  structure.  Clearly,  beyond  this  point 
one  has  to  make  subjective  decisions  regarding  the  probability 
of  the  occurrence  of  maximum  force.  Figures  similar  to  16a 
through  18b  enable  one  to  predict  the  force  coefficients 
which  would  have  prevailed  under  less  favorable  conditions 
(reduced  coherence,  directional  waves,  body  interference, 
etc. )  . 
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VII.  CONCLUSIONS 


The  forces  acting  on  yawed  smooth  and  rough  circular 
cylinders  in  sinusoidally  oscillating  uniform  flow  have  been 
extensively  investigated,  and  the  following  conclusions  have 
been  reached: 

1.  The  independence  principle  does  not  apply  over  signif¬ 
icant  portions  of  the  range  of  K  values  investigated.  The 
drag  and  inertia  coefficients  differ  significantly  from  those 
deduced  from  the  90-degree  case  for  all  45-degree  and  60- 
degree  smooth  and  rough  cylinder  cases. 

2.  The  effect  of  roughness  is  to  bring  the  yawed  and 
unvawed  cylinder  force-coefficients  closer,  except  in  the 
drag-inertia  dominated  regime. 

3.  Yaw  has  a  very  strong  effect  on  the  coherence  of 
vortices  and  hence  on  the  transverse  force  coefficient, 
particularly  in  the  drag-inertia  dominated  regime.  The  lift 
coefficient  is  largest  for  the  90-degree  cylinder,  as  might 
have  been  expected,  and  decreases  sharply  with  yaw. 

4.  An  extensive  sensitivity  analysis  has  shown  that  the 
force-transfer  coefficients  are  affected  by  the  phase  shift 
between  force  and  velocity  to  various  amounts  dependent  on 
the  flow  regime  and  the  magnitude  of  the  phase  shift. 

5.  Finally,  the  results  have  shown  that  the  force 
coefficients  and  the  validity  of  the  Morison  equation  are 


sensitive  to  disturbances  imposed  on  the  in-line  force, 
particularly  in  the  frequency  range  where  the  disturbance 
frequency  matches  the  wave  frequency.  The  stability  envelope 
has  been  delineated  through  extensive  calculations  to  enable 
one  to  interpret  the  scatter  of  the  data  obtained  in  the 
ocean  environment. 
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Figure  2a.  C  Versus  K  for  3  =  4057,  a  =  60  Deg.,  k/D  =  0.00. 
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Figure  2b.  C  Versus  K  for 
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Figure  3a.  C,  Versus  K  for  B  =  4057,  u  =  60  Deg.,  k/D  =  0.00. 
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Figure  3b.  C.  Versus  K  for  B  =  2271,  a  =  60  Deg.,  k/D  =  0.00. 


2.00 


Figure  3c.  C,  Versus  K  for  &  =  1021,  u  =  60  Deg.,  k/D  =  0.00. 
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Figure  3d.  C,  Versus  K  for  6  =  2322,  a  =  60  Deg.,  k/D  =  0.01. 


Figure  3e.  C  Versus  K  for  n  =  1042,  <*  =  60  Deg.,  k/D  =  0.01. 


Figure  4a.  carms  Versus  K  f°r  &  =  4057,  a  =  60  Deg.,  k/D  =  0.00. 
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Fiqure  4b.  C  Versus  K  for  8  =  2271,  a  =  60  Deg.,  k/D  -  0.00. 
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Figure  4e.  C  „  Versus  K  for  6  =  1042,  a  =  60  Deg.,  k/D 
arms 


Figure  5b.  C  Versus  K  for  f3  =  2271,  a  =  60  Deg 


Figure  5c.  C  Versus  K  for  B  =  1021,  a  =  60  Deg.,  k/D  =  0.00 
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Figure  5d.  Cmj_  Versus  K  for  e  =  2322,  a  =  60  Deg.,  k/D 
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Figure  5e.  Cmr  Versus  K  for  3  =  1042,  a  -  60  Deg.,  k/D  0.01 
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Figure  6a.  C,  Versus  K  for  6  =  4057,  u  =  60  Deg.,  k/D 
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Figure  6b.  C  Versus  K  for  t>  =  2271,  u  =  60  Deg 


Figure  6c.  C Versus  K  for  6  =  1021,  u  =  60  Deg.,  k/D  =  0.00. 


Figure  6d.  C,  Versus  K  for  fi  =  2322 
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Figure  6e.  C,  Versus  K  for  f i  =  1042,  a  =  60  Deg.,  k/D  =  0.01 
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Figure  7b.  ^irrns  Versus  K  for  8  =  2272,  a  =  60  Deg.,  k/D 


Figure  7c.  C.  Versus  K  for  £S  -  1022,  a  =  60  Deg 
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Figure  8a.  Comparison  of  C  Versus  K  for  fi  =  4040,  k/D  =  0.00, 
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Fiyure  8b.  Comparison  of  C  Versus  K  for  8  =  2345,  k/D  =  0.00, 


Figure  8c.  Comparison  of  C  Versus  K  for  ft  -  1045,  k/D  -  0.00, 
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Figure  8e.  Comparison  of  C  Versus  K  for  fi  =  2325,  k/D  =  0.01, 
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Figure  9c.  Comparison  of  C-.  Versus  K  for  6  =  1045,  k/D  =  0.00 
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Figure  9e.  Comparison  of  C,  Versus  K  for  3  =  2325,  k/D  =  0.01 
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Figure  10b.  Comparison  of  C  Versus  K  for  6  =  2345,  k/D  =  0.00 

a  rius 

a  =  90,  60,  and  45  Deg. 
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Figure  10c.  Comparison  of  carms  Versus  K  for  3  =  1045,  k/D  =  0.00, 
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Figure  lOd.  Comparison  of  Carms  Versus  K  for  3  =  4100,  k/D  =  0.01 
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Figure  lOf.  Comparison  of  Car;ns  Versus 
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Comparison  of  C  Versus  K  for  3  —  1040,  k./D  -  0.01 


Figure  12c.  Comparison  of  Cdr  Versus 
a  =  90,  60,  and  45  Deg. 
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Figure  12d.  Comparison  of  Cdr  Versus  K  for  M  4100,  k/D  -  0.01, 
=90.  60,  and  45  Deg. 
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Figure  12f.  Comparison  of  C,  Versus  K  for  e  =  1040,  k/D  =  0.01, 
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Figure  13c.  Comparison  of  C^rmg  Versus  K  for  8  =  1060,  k/D  =  0.00, 


Figure  13e.  Comparison  of  C,  Versus  K  for  e  =  2390,  k/D  =  0.01 
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Figure  14a.  Force-Transfer  Coefficients  Versus  Phase  Angle  for  K  -  6.66. 
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Figure  14d.  Force-Transfer  Coefficients  Versus  Phase  Angle  for 
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Figure  15a.  Normalized  Force  Plus  a  Random  Disturbance  of  Frequency 
Ratio  1. 
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Figure  15b.  Normalized  Force  Plus  a  Random  Disturbance  of 


Figure  15d.  Normalized  Force  Plus  a  Random  Disturbance  of  Frequency 
Ratio  36. 
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Figure  15e.  Normalized  Force  Plus  a  Random  Disturbance  of  Frequency 
Ratio  180. 
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18b.  Response  of  Cd  and  the  Corresponding  Cm  to  Random 

Disturbances  of  Various  Frequency  Ratios,  K  =  25.44 
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